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Uniqueness, up to a global phase, is proved for phasing with a random roughly 
known mask (RKM) of high uncertainty. Phasing algorithms alternating 
between the object update and the mask update are systematically tested and 
demonstrated to have the capability of nearly perfect recovering the object 
and the mask (within the object support) simultaneously for up to 50% mask 
uncertainty. The uncertainty threshold for successful recovery can be further 
lifted to more than 70% under the additional sector condition on the object. 
The phasing algorithms with RKM are robust to low resolution masking as 
well as other types of noise. @ 2012 Optical Society of America 

1. Introduction 

Phase retrieval is the problem of reconstructing an unknown image from its Fourier intensity 
data and is fundamental in many applications. Recent breakthroughs center around diffrac- 
tive imaging of non-periodic objects, combining the penetration power of hard X-ray and 
the high sensitivity of lensless imaging |1-3|. 

Despite tremendous progresses many questions, fundamental as well as algorithmic, remain 
to be solved. For example, the standard phasing algorithms, based on alternating projections 
[4]-[6|, are plagued by stagnation and spurious errors partly due to intrinsic non-uniqueness 
of the standard phasing problem. 

We believe that the two problems, non-uniqueness and non-convergence, can be solved 
in one stroke. The present work is a sequel to our previous paper |8|. A centerpiece of our 
approach is the random-mask phasing method J7||8j which is shown to yield a unique solution. 
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up to at worst a global phase factor^ as well as superior numerical performances, including 
rapid convergence, much reduced data and noise stability. In particular, the random mask 
method is robust to various types of noise, including Gaussian, Poisson and mask noises, 
with a noise amplification factor about 2. 

A critique that can be leveled against the random mask approach is the use of exact 
knowledge of the mask which is not always available. Here we present an approach specifically 
to the problem of high mask uncertainty that goes well beyond the stability to mask errors 
demonstrated already in |8|. We will show that nearly perfect recovery can be achieved for 
more than 50% mask uncertainty. 

Our approach is based on the stronger notion of uniqueness (i.e. up to a global phase 
factor)^ and we extend two uniqueness results of |7| to the present setting with roughly 
known mask (RKM). Instead of running phasing algorithms with a fixed erroneous mask, we 
design algorithms to recover the object and the mask simultaneously. At each iteration, the 
object and the mask are updated alternatively aiming at fitting the object constraint, the 
mask constraint as well as the Fourier intensity data. As shown below our numerical schemes 
can accurately recover the object with close to 50% mask errors and, under an additional 
sector condition, more than 70% mask errors. 

The effect of a mask amounts to changing the original object / to the masked object 

(?(n) = A(n)/(n), n e (1) 

where A is an array representing the mask and d is the dimension. In the standard phasing 
problem the uniform mask (UM) with A = 1 is used. In our approach the mask A is random 
and only roughly known. For ease of discussion, we will focus on the case of phase masks 

A(n) = exp (^(/)(n)), (/)(n) G [0, 27r], n G (2) 

We assume that the true phase angles 0(n), in radian, of the unknown mask lie within Sn 
from the given initial estimates (/>o(n) for all n, i.e. 

0(n) G |0o(n)-57r,0o(n) + M- (3) 
Here and below we adopt the following notation: G |a, 6] means 

J a(mod 27r) < 6'(mod 27r) < 6(mod 27r) if a(mod 27r) < 6(mod 27r) 

\ a(mod27r) < 6'(mod27r) < 27r or < 6* < 6(mod27r) else 

For ease of notation, we shall write below the inequality (pi) as 0(n) G |0o(ii) T <^7r]. 
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The paper is organized as follows. We state the uniqueness theorems for phasing with a 
random RKM in Section [2] and give the proofs in Appendices A and B. We discuss the basic 
algorithm of Alternating Error Reduction (AER) and prove the residual reduction property 
in Section |3] and Appendix C. We discuss the Douglas-Rachford-Error- Reduction (DRER) 
algorithm in Section [4] and the algorithms with two sets of Fourier intensity data in Section 
[5) We present numerical results in Section |6] and conclude in Section [7[ 

2. Uniqueness 

Let n = (ni, . . . , rid) ^ and z = (zi, . . . , z^i) G C^. Define the multi-index notation = 
z'^^Z2^ . . • . Let C{N) denote the set of finite complex- valued functions on 1/ vanishing 

outside A/" = {0 < n < N},N = (A^i, A^2, • • • , A^d). Here m < n if < Uj^j. Set 

d 

W\ = \[{N, + l). 

The z-transform F(z) = ^n/(^)^~" of / ^ ^(-^) analytic continuation, from the 

d-dimensional torus, of the the Fourier transform for z = (exp {2Tii(jJi)^ . . . , exp (27ricj(i)), cjj G 
[0,1]. 

If the absolute value of the Fourier transform is sampled on the lattice 



then the autocorrelation function C/(n) = Ylm^M f + ^)/(^) is uniquely determined. 
Sampling on C produces the oversampling ratio a ^ 2^ where 

Fourier intensity data number 
unknown image pixel number 

However, the uniqueness of the autocorrelation given the Fourier intensity data does not 
imply the uniqueness of the object. 

First, there are three types of global ambiguities: (a) constant global phase, /(•) — > 
exp (i6^)/(-), for some 9 G [0,27r], (b) spatial shift, /(•) — > /(• + m), for some m G Z^, 
(c) conjugate inversion, /(•) — > /(N — •). Conjugate inversion produces the so-called twin 
image. 

Our basic tool is the following improved result of irreducibility [9, 10 with, however, a 
more practical and useful perspective. 

Proposition 1. |7| Let f be a finite complex-valued object of rank > 2. Let {A(n)} be 

continuous random variables on nonzero real algebraic varieties {V(n)} in C{c^ M^) with 
an absolutely continuous joint distribution with respect to the standard product measure on 
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rinGS '^(^) ^here T> CN^ is the support set of f. Then the z-transform of the masked object 
is irreducible up to a power of with probability 1. 

The main point here is that while the classical result (9,10 works for generic (thus random) 
objects from a certain ensemble Proposition [T] deals with a given, deterministic object of rank 
> 2. This improvement is achieved by endowing the probability measure on the ensemble of 
masks, which we can manipulate, instead of the space of objects, which we can not control, 
as in the standard setting. 

As a consequence of Proposition [!} the global ambiguities are the sole ambiguities possible 
as far as the masked object ([T]) is concerned p!o|. 

This is not all. We are able to remove all ambiguities, with the only exception of a global 
phase, for the original object / even when the random mask is only roughly known. 

Theorem 1. Let f be a real-valued object of rank > 2. Suppose the exact mask phases {(/)(n)} 
are independently and uniformly distributed on [0,27r]. Suppose the uncertainty of the mask 
estimate Aq = {exp {i(j)o{n))} in ^ is S < 0.5. 

Suppose that another non-negative image f and mask estimate A = {exp (i0(n))} satisfying 

0(n)e[0o(n)T<^7rl (5) 

together produce the same Fourier intensity data on C as do f and A. Then, with probability 
no less than l — \J\f\{25)^^/'^^, /(n) = /(n) Vn and furthermore 0(n) = + (/)(n) for a constant 
^ G M wherever /(n) 7^ 0. Here a positive number x, [xj is the greatest integer at most x. 

Remark 1. // the object is known to be non-negative, then 5 can be any number in [0, 1) 
and uniqueness holds with probability no less than 1 — lA/'l^^'^^^^ . 

For complex- valued objects uniqueness requires two independent sets of Fourier intensity 
data. 

Theorem 2. Let f be a complex-valued object o/rank > 2. Let the first mask A^-^^ = X in 
Theorem^ with the initial mask estimate Aq satisfying 

Suppose the second mask A^^^ is exactly known and the z-transform ofX^^^f is irreducible up 
to a power of z~^ . Moreover, assume the non- degeneracy condition that there is no m ^ Q 
such that A^^^(n + m)/(n + m) = exp exp (ir7(n))A*^^^(n)/(n), Vn, or no m such that 
A(2)(m — n)/(m — n) = exp exp (i77(n))A(^^(n)/(n), Vn, for some ^ G [0, 27r], |ry(n)| < 

Suppose that for a phase mask A with and an object f the two pairs of masked 
objects Xf and Xf, X^^^ f and X^'^^f, respectively, produce the same Fourier magnitudes 
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on C. Then with probability no less than 1 — |A/'|5L'S'/2J jf^^^^ ^ exp (wi)/(n), Vn^ and 
A(n) = exp (iz/2)A(n) if f{n) ^ 0^ where V\^V2 ^ [0,27r]. 



Remark 2. Clearly, most objects and masks obey the non- degeneracy condition. 

The proofs of Theorems [T] and [2] are given in Appendices A and B, respectively. 

Both theorems assert that not only the uniqueness of the object but also the uniqueness 
of the mask, up to a constant phase, inside the object support. This motivates the design of 
our numerical schemes that turn out to be capable of enforcing uniqueness and hence nearly 
perfect recovery even in the presence of high mask uncertainty. 

3. Alternating Error Reduction (AER) 

Let A be the diagonal matrix with diagonal elements {A(n)} and let $ represent the discrete 
Fourier transform. Denote the Fourier intensity data vector by y = |$A/| where A/(n) = 
A(n)/(n). 

3. A. Object Update 

Given the object estimate fk and mask estimate A^^ at the k-th iteration, we use standard 
phasing algorithms to obtain fk+i- 

Let O denote the ensemble of objects / satisfying various object constraints (see Section 
[6]). Let Vo be the orthogonal projection onto O (cf. |8| for details about numerical imple- 
mentation of Vo) and Vi^k = A^^$~^T$A/e, where T is the intensity fitting operator 

Tr(, A-i ^Mexp(^ZG(a;)) if |GM| > 

^^^'^^-{fM if|GM| = o- 

Here and below G [0, 27r) denotes the wrapped phase angle of z. When z = 0, is taken 
to be unless specified otherwise. 

Error Reduction (ER) takes the form f^+i = VoVi^kfk which can be conveniently repre- 
sented by the diagram in Figure [l](a). 

Let r{f^X) = II |$A/| — 1^ || denote the residual. With a phase mask, ER enjoys the 
residual reduction property [sl: 



r{fk+iAk) < r{fkAk)' (7) 
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//C+15 fk 

Vo 



9k 



9k 



(a) object update 



^k+l^ ^k 

Qm 



<i)-ir<i> 



(b) mask update 



f^k 



Fig. 1. Alternating Error Reduction (AER) between object and mask. 



3.B. Mask Update 

Based on the newly updated object estimate fk+i, we define Qi^k as 

r $-ir$A,/fc+i(n)//,+i(n) if /,+i(n) ^ 
a; = Qf,,A,(n) = j . (8) 

[ Afc(n) else 

Let be the ensemble of phase masks satisfying the phase uncertainty constraint ([5]): 

M = {A|Vn,ZA(n)e[(/.o(n)T<^7rl}. (9) 

Let Qm be the orthogonal projection onto Ai. Note that can be computed pixel by pixel 
as follows. 

Let a = (00 (n) — 57r)(mod 27r), b = ((/>o(n) + 57r)(mod 27r) and 



c = 



TT + (a + 6)/2 (mod27r), if a < 6 
(a + b)/2 (mod27r), else. 



Then can be expressed as 



exp(zZAUn)) if ZA^n) G [a, 61 
QmKin) = { exp (lb) if ZA^n) G [6, cj (10) 

exp(m) if ZA'^(n) G |c, a]. 

Since the the object and the mask have interchangeable roles, we set A/e+i = QmQi,k^k 
in the spirit of ER (see Fig. [l](b)). Note the differences between the mask update rule here 
and that of the extended ptychographical engine (ePIE) ((4) in Ref. |11|): First, ([s]) uses the 
newly updated object //e+i while ePIE uses the previous one. Second, more importantly, the 
rough prior knowledge about the mask is enforced by Qm here while ePIE does not consider 
this aspect. 

Now we prove the following residual reduction property. 
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Lemma 1. With we have 

The proof of Lemma [T] is given in Appendix C. 
Define the Alternating Error Reduction (AER) as 

(A+l,A,+i) = {VoPf^kfk.QmQfM^ (11) 

In words, AER alternates between updating the object and the mask estimates until the 
iteration converges. 

Lemma [l] and ([T]) together yield the following residual reduction property for AER. 

Theorem 3. AER ( fTTp has the residual reduction property: r(//e+i, A^^+i) < r{fk^Xk)- 

In our numerical tests, we find that while works well for real- valued objects, for 
complex- valued objects the following alternative rule works better 

n Vfr.^ - / e^P(^^^Un)) if /AUn) G [a, 61 
[ Ao(n), else. 

With Qm we have the alternative version of AER 

(A+l, A/e+l) = {VoVi^kfk, QuiQi^k^k)- (13) 

4. DRER 

AER (either version) by itself converges slowly, typically taking up to several thousands 
steps for accurate recovery in our numerical tests. To speed up convergence we consider the 
Douglas Rachford (DR) algorithm also called the averged successive refiections (ASR), 
which coincides with the hybrid input-output (HIO) algorithm for the parameter /3 = 1 when 
object value constraints are not imposed p!3] 



h» = '-^^^h (14) 



where TZq = 2Vo — /, TZf^k = 27^f,/c — I are refiection operators. 
Define the DRER iteration as 



(A+i,A,+i) = Q(/ + 7^o7^f,,)^,Q^Qf,,A,^ (15) 
and the alternative version as 

(A+i,A,+i) = Q(/ + 7^o7^f,,)^,S^Qf,,A,y (16) 
To strictly enforce the mask constraint, we do not use DR for mask update. 
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5. AER/DRER with two sets of data 

Let A^-^^ = A and A^^^ be two masks with which two sets of Fourier intensity data Y = |$A/| 
and y(^^ = I^A^^^/I are measured on C. Let T and T^^^ be the intensity fitting operators 
corresponding to Y and Y^'^\ respectively. 

For simphcity of presentation we assume the second mask (random or deterministic) is 
exactly known and independent from the first mask which is random. In this case, there is 
no need for the second mask update. 

Suppose fk and Xk are the image and the mask recovered at the end of the k-th iteration. 
At the {k + l)-st iteration, the image is first updated from fk to fk+i based on A^^^ and A^^^. 
Let Vk = {Ak)-^^-^T^Ak and V^^^ = (A(2))-i$-ir(2)$A(2). 

The AER and DRER algorithms with two masks are respectively 

ifk+iAk+i) = {VoV^^^Vkfk,QmQi,kXk), k = o,i,--- (17) 
ih+iAk+i) = Q(/ + 7eo(2p(2)p, _/))/,, Q^Q,,Afc) . (18) 

The following alternative versions work better for complex-valued objects: 

{fk+i,Xk+i) = (VoV^^^Vkfk,QmQi,kXk), k = 0,l,--- (19) 
{fk+i,Xk+i) = Q(/ + 7eo(2p(2)p, _/))/,, Q^Qf .A,^ (20) 

6. Numerical Simulations 

In this section, we test our numerical schemes by performing phasing with a RKM. 

The original images are the 256 x 256 cameraman and the 138 x 184 phantom. We sur- 
round both images by dark (i.e. zero- valued) borders to create images of loose support. The 
resulting images to be recovered are 282 x 282 cameraman and 200 x 200 phantom of loose 
supports (Fig. [2]). Objects of loose support are harder to recover than objects of tight sup- 
port. Everything else being the same, the larger the dark borders the harder it is to recover 
the objects. 

We set the true phases of the RKM to be (/)(n) = (/>o(n) + t(5, n) where t(5, n) are inde- 
pendent, uniform random variables in [—7t6^7t6] with uncertainty S G [0, 1). 

The mask is called a high resolution mask (HRM) if 00 (n) are independently and uniformly 
distributed in [0, 27r] (assumed in Theorems [l] and [2]) and a low resolution mask (LRM) if 
00 (n) are not independently distributed but grouped in independently distributed 10 x 10- 
blocks, each of which is a (random) constant. For our examples, a LRM has less than 900 
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(a) 282 X 282 cameraman 



(b) 200 X 200 phantom 



Fig. 2. Test images of loose support. 



and 400 degrees of freedom, respectively, for the cameraman and the phantom. A HRM, 
however, always performs better than a LRM as HRM has many more (100 times) degrees 
of freedom (cf. Fig. g). For 200 x 200 HRM and 20 x 20 LRM used for Fig. g see Fig. g 
where the gray level represents the mask phases in [0, 27r]. 

When a second set of Fourier data is used (for complex- valued objects), it is generated 
withaUM (i.e. A^^) = 1). 

A standard way to utilize the oversampled data for a > 1 is to enlarge the original image 
by adding corresponding number of zero pixels which is then enforced as an additional object 
constraint. This procedure is called the oversampling method and implemented in all 
our simulations with a = 4. 

6. A. Error and Residual 

To estimate the recovery, we define relative error and relative residual as follows. Let / and A 
be the recovered image and mask respectively. The relative error of reconstruction is defined 
as 



e{f) 



if absolute uniqueness holds 
min 11/ — exp (^^)/||/||/|| if uniqueness holds only up to a global phase 



Let A be the diagonal matrix whose diagonal elements are A(n). The relative residual is 
defined as 

, ^ II / - |$APo/| II 
P[f, A) = ^ 



(a) 200 X 200 HRM 



(b) 20 X 20 LRM 



Fig. 3. High (a) and low (b) resolution masks 



where Vo is introduced to enforce the object constraints in the case of DRER. 



6.B. Performance study of AER and DRER 



First we use AER (11) to recover the non-negative images with the stopping rule \\fk+i — 
fk\\/\\fk\\ < 0.05% and one LRM of uncertainty S = 0.2. The results, shown in Fig. [4| are 
noisy and inaccurate with 33.75% error for the cameraman and 44.74% error for the phantom. 
Consistent with the residual reduction property (Theorem |3]), the residual curves in Fig. [4] 
are monotonically decreasing. 

Much improvement can be gained by running DRER, followed by AER. For real-valued 

fk 
fk 



objects, we use the version of DRER (15). DRER (15) is stopped when \\fk+i 



1%, with the maximum of 400 steps, and AER (11) is terminated when \\fk+i 
0.05%, with the maximum of 400 steps. As shown in Figure [5| for the cameraman, 131 DRER 
and 6 AER steps took place producing 1.43% error while for the phantom, 75 DRER and 5 
AER steps took place with 0.33% error. Consistent with Theorem [!} the mask errors occur 
only outside the object supports. 

Next we consider the case of the generic complex- valued objects with one UM and one 



LRM of uncertainty S = 0.2. We apply the alternative versions of DRER (20) and AER (19) 



which tend to produce better results than ( |15[ ) and ( |11[ ) for complex- valued objects. DRER 

-fk\ 



(20) is stopped when H/^+i 



fk\\ < 1%, with the maximum of 400 steps, and AER 
( [19| ) is terminated when H/^+i — All/ 1| All < 0.05%, with the maximum of 400 steps. Fig. 
[g] shows the results for object phases randomly distributed on [0, 27r]. Both algorithms ran 
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their full course of 400 steps with 2.39% error for the cameraman and 1.37% error for the 
phantom. The mask errors occur only outside the object supports, consistent with Theorem 

El 

Still next, we let the object phases be randomly distributed on [0,7r/2] (the sector con- 
straint) and include this prior knowledge into Vo- With the additional sector constraint, 
we found that the following stopping rule can significantly reduce the number of iterations: 



DRER (20) is stopped if the residual increases in five consecutive steps, with the maximum of 
400 steps, and AER ( [T9| ) is terminated when H/^+i — /fc||/||//c|| < 0.05%, with the maximum 
of 400 steps. Fig. [7| shows the results with one UM and one LRM of uncertainty 5 = 0.2. 
With the new stopping rule and the sector constraint, 24 DRER and 400 AER steps took 
place with 1.95% error for the cameraman and 24 DRER and 365 AER steps with 0.61 % 
error for the phantom. 

Finally Fig. [s] shows the averaged relative error e(/), after 5 runs of independently chosen 
initial guesses for the object, with or without mask update, as a function of the mask uncer- 
tainty of HRM or LRM for (a) non-negative phantom, (b) generic complex- valued phantom 
and (c) phantom under the sector condition. We use the same stopping rules and updat- 
ing rules as above for each case, except that the maximum number of steps is changed to 
200 + 5 • 1000 for DRER and AER separately to deal with variable uncertainty. 

Without mask update the error curves are roughly linear with the noise amplification 
factor roughly 2 (blue and black curves), consistent with our previous results reported in |8|. 
With mask update, the results (pink and red curves) are drastically improved for uncertainty 
up to 50% (Fig. [s] (a) and (b)) and, under the sector constraint, nearly 70% uncertainty in 
the mask phases (Fig. [S] (c)). 

7. Conclusion 

We have proved the uniqueness, up to a global phase, for phasing with RKM with probability 
exponentially close to one, depending on the object sparsity and mask uncertainty. We 
have designed algorithms that achieve nearly perfect recovery for mask uncertainty up to 
half of that promised by the uniqueness results. Under the additional sector condition, the 
algorithms can achieve nearly perfect recovery for up to 70% uncertainty in the mask. As 
a by-product of object recovery the unknown mask can be recovered accurately within the 
object support. The numerical performance is robust to low resolution mask as well as other 
types of noises as shown in the previous study [sj. 

Our method can be easily extended to general masks as long as they modulate object 
phases randomly (i.e. they are not just amplitude masks). This can be done by first re- 
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constructing not the original object but the object multiphed with the intensities (assumed 
nonzero) of the mask by using the normahzed (phase) mask. In the second step, the original 
object is recovered by dividing by the intensities of the mask. 

A. Proof of Theorem [T] 

Proof. In view of Proposition [l| there exist some m and ^ G [0, 27r] such that either 

A(n)/(n) = exp (i^)A(m + n)/(m + n) (21) 



or 



A(n)/(n) = exp {i9)X{in — n)/(m — n). 



(22) 



In the case of (21) with any m 7^ and any 9 G [0, 27r], 

. ...JA(m + n)|exp(z(/)(n + m)) 

/(n) = exp (lO) /(n + m). 

|A(n)| exp (zZA(n)) 

Consider the lS/2\ independently distributed r.v.s. of A(n + m) where /(n + m) 7^ corre- 
sponding to [5^/2] nonoverlapping pairs of points {n, n+m}. For every n where /(n+m) 7^ 0, 
a proper choice of ZA(n) makes /(n) real- valued if and only if either 



cPin + m) e IM^) - - If in + m) t ^ttJ 



or 

(/)(n + m) G [(0o(n) - ^ - Z/(n + m) + tt) T ^tt] 

However, (/)(n + m) is independently and uniformly distributed in [0, 27r], so it falls in these 
two regions with probability 25. The probability for every such /(n) to be real- valued is 25 
and hence the probability for all /(n) with m 7^ to be real- valued is at most (25)^*^/^^. 

The union over m 7^ of these events has probability at most \J\f\{25)^^/'^^ . Therefore, 
with probability at least 1 — \J\f\{25)^^/^^ , m = and exp (i0)A(n)/(n) = A(n)/(n) Vn which 
further implies that /(n) = ±/(n) Vn and A(n) = ±exp (i0)A(n) on n where /(n) 7^ 0. 



Likewise the probability for all /(n) given by (22) to be real- valued for any m is at most 



|A^|(25)L5/2J. □ 
B. Proof of Theorem [2] 

Proof. In view of Proposition [l} for some mi, m2 and 6*1, ^2 ^ [0, 27r] either 

exp (^^l)A(n + mi)/(n + m^ = A(n)/(n) (23) 
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or 



as well as 



or 



exp (i6'i)A(mi - n)/(mi - n) = A(n)/(n) 
exp (i^2)A(')(n + m2)/(n + 1112) = X^^\n)f{n) 



exp (i^2)A(2)(m2 - n)/(m2 - n) = X^^\n)f{n). 



(24) 
(25) 
(26) 



There are four possible combinations of (23), rt24L i\25h and (26) 



In the case of rt23P&(25), we have 



exp (i^i)A(n + mi)A(2)(n)/(n + mi) = exp {i92)X^^\n + m2)A(n)/(n + ms). (27) 

For any mi ^ and any ^i, ^2 ^ [0, 27r], consider the [5^/2] pairs of independently distrib- 
uted r.v.s. of A(n + mi) where /(n + mi) 7^ corresponding to [5^/2] non overlapping sets 
of points {n, n + mi}. For every n, a proper choice of A(n) makes ([27|) true if and only if 



(/)(n + mi) 

G [((/)o(n) + 02-91 + (^^'\n + m^) - (l)^'\n) + Z/(n + m^) - Z/(n + m^) + SnpS) 

where 0(2) (n) = lX^^\n). 

Since 0(n + mi) are independently and uniformly distributed in [0, 27r], ([28]) holds for each 



n with probability at most S and hence (27) holds for all n at once with probability at most 

The union over mi 7^ of these events has probability at most |A/'|5L'S'/2J ^ Therefore, with 
probability at least 1 — |A/'|5L'S'/2J ^ ^ q ([27l) becomes 



A(n) A(2)(n + m2)/(n + m2) 

= exp {i02 — lOi) — 



A(n) 



A(2)(n)/(n) 



Moreover, if A2/ satisfies the non-degeneracy condition, then m2 = 0, /(n) = 
exp (i6^2)/(n), Vn, and A(n) = exp (i6^i — i6^2)A(n), if /(n) ^ 0, with probability at least 
1-|A/'|5L^/2J. 



In the case of (p3|&(26), we have 



(29) 



exp {i9i)\{xi + mi)A(2^(n)/(n + mi) = exp (i02)A(n)A(2)(m2 — n)/(m2 — n). 

The same argument applies and mi = with probability at least 1 — |A/'|(^L'S'/2J ^ ( [2^ 
becomes 

A(n 



A(n) 



exp {162 — iOi) 



A(2)(m2 - n)/(m2 - n) 
A(2)(n)/(n) ■ 
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which violates the non-degeneracy condition. In other words, (23)&(26) holds with proba- 
bility at most \Af\S^^/^K 



□ 



Similar conclusions follow in the case of (24)&(25) and (24)&(26). 
C. Proof of Lemma [T] 

Proof. Since the operator T enforces the measured Fourier intensities Y 

r{fk+i,Xk+i) = II |$Afc+i/fc+i| - F II 

< ||$Afc+iA+i - r$Afc/fc+i|| = ||Afc+iA+i - $-'r$Afc/fc+i|| 

by the unitarity of the Fourier transform. By splitting the summation and using the definition 
(|8|, the rightmost term becomes 

1/2 



J2 |A+ir(n) A,+i(n)-AUn) 

v/fc+l(n)7^0 



/fe+l(n)=0 



. (30) 



Now since Ajt+i(n) = Q^yj^{n) is a pixel- wise projection of A'^(n), \Xk^i{n) — A'^(n)| < 
|Ajt(n) — A'^(n)| and hence ( [30| ) is less than or equal to 

1/2 

< I IA+ir(n) A,(n)-AUn)% J] \^-'T^A,h^,{n) 

\/fe+i(n)/0 /fc+i(n)=0 

= \\^kfk+i - *"^T$AfcA+i|| = \\^Akfk+i - T^^kfk+i\\ 
= II \^Akfk+i\ -y \\= r{fk+u Afc) 



which is the desired result. 



□ 
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(a) error = 33.75% 



(b) 



(c) res = 6.33% 




Fig. 4. Recovery of non-negative images by AER with one LRM of 5 = 0.2. 
The middle column shows the absolute phase differences between A and A. The 
right column shows the relative residual at each iteration. 
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(a) (b) (c) 




Fig. 5. Recovery of non-negative images with one LRM of 5 = 0.2. (a) the 
recovered cameraman / by 131 DRER + 6 AER steps. e(/) ^ 1.43% and 
A) ^ 0.25%. (d) the recovered phantom / by 75 DRER + 5 AER steps. 
e(/) ^ 0.33% and p{f^X) ^ 0.12%. The middle column shows the absolute 
phase differences between A and A. The right column shows the relative residual 
at each iteration. 
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(a) 



(b) 



(c) 




(d) 



(e) 



re residual at each iteration 



(f) 



Fig. 6. Recovery of generic complex-valued images with one UM and one LRM 
of 6 = 0.2. (a) absolute values of the recovered cameraman / by 400 DRER 
+ 400 AER steps. e(/) ^ 2.39% and A) ^ 1.03%. (d) absolute values 
of the recovered phantom / by 400 DRER + 400 AER steps. e(/) ^ 1.37% 
and A) ^ 0.65%. The middle column shows the absolute phase differences 
between A and A. The right column shows the relative residual at each iteration. 
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(a) 



(b) 



(c) 





50 100 150 200 250 300 350 



(d) 



(e) 



(f) 



Fig. 7. Recovery of complex- valued objects under the sector constraint with one 
UM and one LRM oi S = 0.2. (a) absolute values of the recovered cameraman 
/ by 24 DRER + 400 AER steps. e(/) ^ 1.95% and A) ^ 0.94%. (d) 
absolute values of the recovered phantom / by 24 DRER + 365 AER steps. 
e(/) ^ 0.61% and p{f/X) ^ 0.37%. The middle column shows the absolute 
phase differences between A and A. The right column shows the relative residual 
at each iteration. 
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(a) (b) (c) 



Fig. 8. Averaged relative error e(/) of 5 independent runs versus the percent- 
age of mask uncertainty for (a) non-negative phantom (b) generic complex 
phantom and (c) phantom under the sector condition. The stopping rules and 
mask updating rules are the same as described in the main text for each case 
with the maximum of 200 + 1000 • S steps for DRER and AER separately. 
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